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Fine and hyperfine interaction on the light cone 
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Abstract. The formalism for the spin interactions in the front form (light-cone) is re-phrased in terms of 
an instant form formalism. It is shown how to unitarily transform the Brodsky-Lepage spinors to Bj0rken- 
Drell spinors and to re-phrase the so called spinor matrix in terms of the interactions one is familiar with 
from atomic and Dirac theory. — One retrieves the (relativistic) kinetic correction, the hyperfine and 
the Darwin term which acts even when wave function is spherically symmetric. One also retrieves angular 
momentum dependent terms like the spin-orbit interaction in a relativistically correct way; and one obtains 
additional terms which thus far have not been reported particularly various L 2 -dependent terms. Since 
the approach includes the full retardation, one gets additional, thus far unknown terms. The differ from 
atomic and Dirac theory, since there only that part of the vector potential is usually included which is 
generated by the atomic nucleus. Quite on purpose, the paper is kept formal. — 

PACS. ll.10.Ef - 12.38.Aw - 12.38.Lg - 12.39.-x 



1 The light-cone integral equation 

This paper number 3 in a row of 3 [HE] on the bound 
state problem in gauge theory P] deals with the technical 
question of how to formulate the fine and hyperfine inter- 
action in the one-body integral 'master' equation which 
has been previously derived |2B3]- 

I therefore jump immediately to Eq.(16) of [2], 



M 2 4> hlh2 (x,'k 1 



my, 



-21 



1 

4^2 



E 



X 1 — X 



i>hih 2 (x,kL) 



dx'(PV ± il> h , h ,(x , ,U ± ) a c (Q) 
= — 



y/x{l-x)x'{l-X') Q 2 

x [ufa, fcijyMfcUi)] [v(k 2 ,h' 2 )^v(k 2 ,h 2 )} . (1) 

Here, M 2 is the eigenvalue of the invariant-mass squared. 
The associated eigenfunction iphih 2 ( x i Hi) is the probabil- 
ity amplitude (x, kj_, hi] 1 — x, — kj_, h 2 \\P q q) for finding the 
quark with momentum fraction x, transversal momentum 
kj_ and helicity hi, and correspondingly the anti-quark. 
Their (effective) masses are denoted by mi and m 2 , and 
u(ki,hi) and v(k 2 ,h 2 ) are their Dirac spinors in Lepage 
Brodsky convention, as given in py. The (effective) cou- 
pling function a c (Q) = |a(Q) is also given in [3J- The 
kernel is governed bv the mean four-momentum transfer, 
Q 2 = \{Q 2 q + Q 2 q ), where 



Q 2 q = -(h - k[) 2 and Q 2 = -(k 2 -k' 2 ) 



(2) 



are the Feynman four-momentum transfers of quark and 
anti-quark, respectively. The regulator function R{Q), fi- 
nally, removes the ultraviolet singularities and regulates 



the interaction. Note that the equation is fully relativistic 
and covariant. It coincides literally with Eq.(4.101) of 3 . 

Krautgartner et al 4 and Trittmann et al [S] have 
shown how to solve such an equation numerically with 
high precision. But since the numerical effort is consider- 
able, it is reasonable to work first with simpler models. 
The aim of the present work is to derive such ones. 

The aspects of regularization and renormalization have 
been emphasized in resulting in an explicit construc- 
tion of the regulator function R(Q). The case was worked 
out within the so called Singlet-Triplet model. Here, I ad- 
dress to go beyond that, particularly to derive a model for 
the spin-orbit interaction, which had been suppressed on 
purpose in [5]. 



2 Transforming the integral equation 

The light-cone integral equation (Tj@) has the unpleasant 
aspect that the integration variables have a completely 
different support, 



< x < 1 , 



< ki < 



Therefore, practically in all of the numerical work partic- 
ularly in rj] and the variable transform 



x(k z ) = 



Ei+k z 



(3) 



k 2 



A 2 , 



with Ei_2 = Ei^(k) — y in 1 2 

has been used to transform to integration variables 

-co < k z < +oo , — oo < kj_ < +oo , 
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with the same support. While k z varies from — oo to +00, 
the x(k z ) varies from to 1. The particles are then de- 
scribed by their front form four-momenta 



k z +E 1 (k), k+ = - k z + E 2 (k), 



k± 2 = - kx , 



k{ =-k z +E 1 (k), k 2 = k z +E 2 (k). 

Or, they are described by the instant form four-momenta 

fc° = E x {k), k° 2 = E 2 (k), 
ki = k , k 2 = - k , 

with k = (kj_,fc 2 ). Such a switching between front form 
and instant form parameterization is possible, since the 
four- vectors of the constituents refer to free particles. The 
free invariant mass of the two particles, 



mf + kf; m\ + k 2 



1 -x 



= (E^k) + E 2 {k)) 2 , (4) 



looks like in the rest frame of the instant form (P = 0). 
For vanishing k it is (mi + m2) 2 . The discrepancy can be 
calculated exactly as |S] 

(E 1 + E 2 f - (mi + m 2 ) 2 = 

(-Ei + E 2 -m 1 - m 2 ) (E 1 + E 2 + mi + m 2 ) , 

and therefore as 

{E 1 +E 2 f - (mi +m 2 ) 2 = 

E? - m? El-mV\ , 

* 1 + 2 2 (Ei +E 2 + mi + m 2 , 
iil + mi b, 2 + m 2 / 

With the reduced mass, 



111 

m r mi m 2 



m\m 2 



mi + m 2 

and the dimensionless A(fe) and B(k), 



(5) 



A(fc) = 



= gx(fc)+g 2 (fc) 

^ E 1 {k)E 2 (k) 



B(k) 



_ -Ei(fc)+m 1 +g2(fc)+m 2 
(1111+1112) 



+ 



(fe)+mi 1 E 2 (k)+m 2 J ' 



(6) 



the free invariant mass and the exact Jacobian of the 
transformation J2J| is therefore 



M f 2 ee = (mi + m 2 y + (mi + m 2 ) — B(fc) , (7) 



dx = x(l — x) 



dk z 1 



A(k) 



(8) 



respectively, see also jH]- 

The transformation from (x, kj_) to (k Z) kj_) will be pre- 
sented in two steps. In the first step, the variables are 



transformed and Eq.Q becomes 



M 



{E 1 {k) + E 2 {k)f ^ hlh2 {k z M 



E 



dk'd 2l k 



, ^x<{\-x')^h> lK (k'M 



y/x(l - X) 



a c (Q) i?(Q) 



Q 2 4ir 2 m r 

One notes that the kernel is not symmetric under the ex- 
change of primed (') and unprimed quantities, as opposed 
to Eq.Q which is symmetric. This asymmetry can how- 
ever be removed as usual, by multiplying the equation in 
a second step with ^Jx(l — x)j y/A(k), 

4>h 1 h 2 (fea.kjj 



m 2 - (Si + E 2 y 



-1ph l h 2 (kz,kL.) = 



E jdk' z d 2 K 

h[,h' 2 

a c (Q) R(Q) 



, VAT 



1 



Q 2 47r 2 m r y/A(k)A(k') 

x mkuhj-fuKhl)] [v(k' 2 ,h' 2 ) lfl v(k 2 ,h 2 )} . (9) 

In the numerical work 4,5 , the reduced wave function 
<j>hih 2 {kz, kj_) is calculated first. It is then converted to 
-ip hl h 2 (x,k±) by 



tl) hl h 2 (x,kj_) = ^ -cj) hlh2 (fc z (x), kj_) , (10) 

v /x(l - x) 

i.e. by the substitution k z — k z (x) which is inverse to 
Eq.©. Krautgartner 0] and particularly Trittmann [5] 
have presented beautiful three-dimensional plots of iphih 2 - 

The variable transformation © is applied here in a 
strict mathematical sense. It does not change the physical 
content, particularly not the eigenvalue spectrum M 2 . The 
transformed integral equation @ looks like an equation 
in usual momentum space. But one should emphasize that 
it continues to be a front form equation with the sole pur- 
pose to generate t/Jhxh 2 (x, kj_). The reduced wave function 
fyhxhi {%■> ki_) has no physical interpretation. 

The kernel of Eq.@ depends on |k'| 2 and Lorenz- 
invariants like Q 2 or [7^] [7^]. It is therefore invariant un- 
der spatial rotations. The occurrence of spin-degenerate 
multiplets in the numerical solutions of |Sj become thus 
understandable e posteriori. It seems as if all light-cone 
specific troubles with rotations of the coordinate system 
[3| are absorbed in the x(l — x)-factor in Ea. (|10|) : The re- 
duced wave function cj)h 1 h 2 (k z , kj_) transforms covariantly 
under rotations while iphih 2 {x,k±) does not. 



3 The Melosh rotated integral equation 

Much of the difficulty in getting the reduced wave function 
in practice [4*||^]. and to understand the structure of the 
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numerical solutions, resides in the so called spinor factor 
(h 1 ,h 2 \S\h' 1 ,h' 2 ) = 

= [^k 1 ,h 1 ) 7 '"u(k' 1 ,h[)}[lj(k 2 ,h' 2 ) ltl v(k2,h2)} . (11) 

The 4x4 matrix in the helicities can be simplified by 

[v(k 2 , h' 2 )j fl v(k 2 , h 2 )} = [u(k 2 , h 2 )^u(k' 2 , h' 2 )} . (12) 

This general spinor identity allows to work only the u- 
spinors. But even then, the spinor factor 



with k r = k x + iky and ki = k x — iky. The Bj0rken-Drell 
spinors are [S]: 



u BU (k,s) 



y/2m(E + to) 



E + to 




o 1 

E + m 

h 

-k z ) 



with E = E(k) as in Eq.(|3J). The four overlap matrix 
elements of Eq. (|Td|) are then calculated as 



(/H,/1 2 |SK,4) = (h u ht\S(x,k ± ;a/ t li_)\h[X) = 

[u(fci,/ii)7 M u(fci,/i' 1 )] [u(k2,h 2 )j u ,u(k 2 ,h 2 )] > ( 13 ) 



(s\w\h) 



y/2p+(E + to) \ k r 



ui —ki 

k+ 



(17) 



m —k\ 
k + + m 



is a terribly complicated function of helicities and light- 
cone momenta x,kj_ and x',k^, as seen in the explicit 

tables in jSj. They become even more complicated if one ^^(s\u;\h) (h\uj^ \s') = 
expresses them in terms of k z , kj_ and k' z , , see [5] . h 

One conjectures that the spinor function 1|13JI is much 
simpler if the Lepage-Brodsky spinors u(k, h) = u LB (k, h) 
are replaced by the Bj0rken-Drell spinors u BD (k, s) 0. As 
to be seen, this is the case indeed. 

The way this can be done was shown first by Krass- 
nigg et al. \7\, and their work shall be repeated here in 
short. - Both, the Lepage-Brodsky and the Bj0rken-Drell 
spinors are solutions to the same equation, the free Dirac 
equation — m) u(p, A) = 0. Hence, they must be linear 
superpositions of each other. We define 



with the rows labeled by s and the columns by h. One 
calculates 



1 



2k+{E + m) 



k + + m ki 
—k r k + + m 



1 
1 



since (p + + m) 2 + pip r = 2p + {E + to), verifying this way 
that the transformation is unitary, indeed. 

Since the spinors appear in bilinear combinations, it is 
convenient to define the unitary direct product 



.LB 



(k,h)=J2u BD (k,s){s\cj\h). 



(sis 2 |ft|/iift 2 ) = (siMM ® (S2M/12) ■ (18) 
Introducing a second reduced wave function (p slS2 by 



(14) 



The transformation matrix (s|w|/i) is independent of the 
Dirac indices a. If both spinors have the same normaliza- 
tion, the transformation is unitary — and then called a 
Melosh rotation PTHIITT] . 

The Lepage-Brodsky spinors are conventionally nor- 
malized as u LB (k, h)u LB (k, h!) = 2mShh', while the Bj0r- 
kcn-Drcll spinors are normalized to unity. We change con- 
vention of the Lepage-Brodsky spinors by requiring 



if SlS2 (k) = J2 (siS2\n\hih 2 )<p hlh2 (k) , (19) 

hi,h 2 

Eq.© can be transformed unitarily to 



M 



(E 1 (k) + E 2 {k)f VsiS2 (k) = 



(mi 



m 2 , 



J2jdk' z d\ 



, <Ps[s' a (k') 



y/A{k)A{k>) < 
x [u(fc 1 ,s 1 ) 7 ' i u(fci,si)] [u(fc 2 ,s 2 h M u(fc 2 ,s 2 )] 



(20) 



u LB (k,h)u LB (k,ti) =8 hh , 



(15) 



Eventually, this implies to multiply the kernel of Eq.© 
with the factor 4toito 2 . 

The expansion coefficients from Ea. l|14(l are then de- 
termined by 



Once one has the wave functions <p Sl s 2 (k z , kj_), one can 
unitarily transform them back to the light-cone wave func- 
tions by 



^A(k z (x),k ± ) 



( S \u J \h)=J2^ D (k,s) u^ B (k,h). 



(16) 



y/x(l - x) 

S2)<p SlS2 {k z (x),kj.) 



(21) 



The Lepage-Brodsky spinors are 0: 



u LB {k,h) = 



1 



k + + to 
k + — to 



-h 

k + + to 
ki 

-k + + m J 



in analogy to Ea. i|10fl . 

The spinors in Ea. (|20|l are Bj0rken-Drell spinors. On 
the technical level, they are much more transparent than 
those of Lepage Brodsky, as to be seen next. The previ- 
ous numerical work done with Ea. (|2()fl would have 
been much easier than with Eq.|0. But at that time, the 
present physical insight was lacking. 
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3.1 Definition of the spinor factor 

The spin dependence in Eq. 12U|) resides in the spinor factor 
(s 1 ,s 2 \S\s' 1 ,s' 2 ) = [S(fci,si)7^it(fci,si)] 



With 7"(1)7 /1 (2) = 7°(1)7°(2) - 7 i (l) 7 i (2), where 
it can be evaluated in closed form with [5j 

X S1 



cr 
-a 



u(ki,si) 
u(k 2 , s 2 ) 



Ei + mi 
2mi 



•k 



£i (k) + mi 



S 2 



m2 



2m 2 



• k 



E 2 {k)+m 2 



has been transformed to the 4x4 matrix of the effec- 
tive light-cone Hamiltonian in usual momentum represen- 
tation, Ea. lf2T)jl . 

/ ( i 3 k'< Sl , S2 |iT eL c(k;k , )|si,4)^(k / ) = AfV SlS2 (k). 

Here and below the explicit summations over the helicities 
are replaced by the Einstein convention. The spectrum of 
invariant mass-squared eigenvalues M 2 is unchanged by 
the transformation. One can subtract a c-number from 
an operator, and divide by a scale, and thus define a 
new Hamiltonian H with new eigenvalues E but the same 
eigenf unctions <p Sl s 2 '- 



(22) #cLc(k; k') = (m a + m 2 f + 2 (mi + m 2 ) iT(k; k') , 

M 2 — (mi + m 2 ) 2 + 2 (mi + m 2 ) E . (24) 



The Pauli spinors are \] = i^j an< ^ ^-l = ( 1^ ' 
The components of the four current = (p, j) are then 



The eigenvalue E has the dimension of an energy and is 
not to be confused with the single particle energy E(k): 



„ -\7n(h c cMl - / E 1 (k)+ mi E 1 (k')+ mi d 3 k' (si , s 2 \ H (k; k') | si , s 2 W s > (k') = Eip Sl S2 (k) , 

Pi — |W(Ki,SiJ7 U(Ki,S 1 )l - y 2^ 2^ J 



X (Si 



1 



(El(fc)+mi) (Bi(fe')+mi) 



(<Ti-k)<ri , cri(o-i-k') 



X (Si 



(Bi(fc)+mi) 1 (£i(k')+mi) 



where H (k; k') = T(k; k') + C/(k; k') . (25) 
The kernels for kinetic and potential energy are given by 
A- 2 



T(k;k') 



2m, 



The spinor factor becomes then most straightforwardly: 
( Sl ,s 2 \S\ S [,s' 2 ) = J EA 2 k l +mi E2{ 2 k l + ; n I (23) 



6^(k-k')5 sl< d S2< B(k), (26) 
1 



(27) 



E 1 (k / )+m 1 E 2 {k')+ni2 
2mi 2 m 2 



Sl, s 2 



(ffl-k)(<Tl-k') 



X 1 + 



(E 1 (k)+m 1 )(E 1 (k')+m 1 ) 
(o-2-k)(er2-k') 



/ (<T1 

V(-EiC 



(B 2 (/c)+m 2 )(B 2 (fc')+m2) 



(fc)+mi) 
(o-2-k)o-2 



cri(g-i-k') \ 

(Ei(k')+ mi ) J 

cr 2 (o-2-k') 



N /A(fc)A(fc' 

respectively. Both A(fc) and -B(fc) were defined in Eq.©. 

5 Explicit calculation of the spinor factor 

This section is devoted to carry out explicitly the mul- 
tiplications in Ea. Q23f) . The bilinear expressions with the 
same cr's can be simplified by means of the identities 



(E 2 (k)+m 2 ) (E 2 (k')+m 2 ) 



s li s 2 



The expression is very much simpler than the long tables 
for the Lepage-Brodsky spinors [B] , indeed. The first two 
lines in the square bracket correspond to the product of 
the charges, pip 2 , and the next two lines to the scalar 
product of the currents, jij2- Note that the plus sign of 
jij2, as opposed to the minus sign in 7 M 7 M = 7°7° — 77. 
This is due to ki = +k and k 2 = — k, see also EQ. (j22f) . 



4 Conversion to a conventional Hamiltonian 

Thus far, the 4x4 matrix of the effective light-cone Hamil- 
tonian in light-cone momentum representation, Eq. (|9"|> . 

J dk' z d 2 ^ (hi , h 2 \H cLC (x, k x ; a/, ki)K , h' 2 ) 



(k-er)cr= k + i(crAk), 
<r(k'-<r)= k'-i(o-Ak'), 
(CT-k)(<T-k 7 ) = kk' + (kAk'). 

One gets thus identically in a first step: 



(Si, 52151s;, S 2 ) - • / M±ga E 2 (k)+m 2 



£i(fe')+mi E 2 {k')+rn 2 
2mi 2m 2 



(28) 



(29) 



91, 32 



X 1 + 



kk ' +iai -kAk / 

(_E 1 (fc)+m 1 )(_E 1 (fc')+m 1 

k-k'+ia-2-kAk' 



(E 2 (k)+ m2 )(E 2 (k>)+77i 2 ) 

_r_ / k+iq-jAk 1 k'— itriAk' \ 
+ ^(EiW+mi) + (Ei(t')+mi)J 

( k+iq- 2 Ak 1 k' — ia 2 f\V. ! \ 
(E 2 (k)+m 2 ) (E 2 (k')+m 2 ) J 
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Carrying out the multiplications one gets in a second step: 
(s 1 ,s 2 \S\s[,s' 2 ) = ^^ E ^ m2 



(30) 



E 1 (k')+m 1 E 2 (k')+m 2 
2mi 2m2 



Sl,S 2 



i k-k , _+ig'i kAk ' , k-k ' +ia 2 -kAk [ 

+ (Bi(fe)+m 1 )(B 1 (fe')+m 1 ) {E 2 (fe)+m 2 ) (E 2 (fe')+m 2 ) 



k-k ' +i<Ji -kAk ' 



k-k'+m2-kAk' 



(E 1 {k)+m 1 )(E 1 (k>)+m 1 ) (E 2 (k)+m 2 )(E 2 (k')+m 2 ) 
k+itJi Ak k+ia'2 Ak i k ' — io\ Ak ' k — ia 2 A]£ 



(Bi(fe)+mi) (B 2 (fc)+m 2 ) ^ (Ei(t')+mi) (B 2 (i')+m 2 ) 

k+itri Ak k 7 — zct 2 Ak / 
(£i(fe)+mi) ' (B 2 (fc')+m 2 ) 



k ' — itTi Ak 7 



k+'to- 2 Ak 



(£?i(fc')+mi) (£ 2 (fc)+m 2 ) 

With Eas. (|32l33l34[l below, one gets in a third step 



( S1 , S2 \s\ S [, s >)=j ^i + ^ g tir 2 



(31) 



2mi 2m-2 



Si, 5 2 
k-k / +i<?2-kAk / 



i + 



+ 



k-k ' -kAk 7 , 

(E 1 (k)+m 1 )(E 1 (k')+m 1 ) "•" (£ 2 0)+m 2 )(,E 2 (fc')+m 2 ) 

k 2 k' 2 -(kAk') 2 -(cri-kAk')(o- 2 -kAk')+i(cri+cr 2 )-kAk'(k-k') 

(£l(fc)+mi)(_Bi(fc')+mi)(-B2(fc)+m 2 )(_E 2 (fe')+m 2 ) 

k 2 -(tT 1 Ak)-(cr 2 Ak) . k' 2 -(o-iAk')-(q- 2 Ak') 



(_Ei(fe)+mi)(_E 2 (fc)+m 2 ) 1 (Ei (k')+mi ) (E 2 (fc' ) +m 2 ] 
kk' + (o-iAk)-((T 2 Ak')+i(<T 1 +(T 2 )-(kAk') 
H (-Ei(fe)+m 1 )(_E 2 (fc')+m 2 ) 

kk ' + (cti Ak ' ) ■ (er 2 Ak)+i(cri +cr 2 ) ■ (k' Ak) 
^ (-Bi(fc')+™i)(-B2(fc)+m 2 ) 

Here, the familiar vector identities like 



a A b • c = a ■ b A c , 

(aAb)Ac = (ac)b -(b-c)a, (32) 
(a A b) • (c A d) = (a • c)b • d (b • c)a • d , 



were used to derive particularly 

(k-k'f=k 2 k' 2 -(kAk'f , 
(k + ioi A k) ■ (k' - ia 2 A k') = 



(33) 



kk 



(kAk') . 



- ia\ A k) ■ (k' - io 2 A k'J = (34) 
+ (<ti A k) • (<j2 Ak') + i {a x + a 2 ) • (i 

In deriving Ea. (|31(l . Ea. l|33|) was used in line 3 of the 
square bracket of the equation. In lines 5 and 6, Ea. (|34|l 
was applied. In line 4, Ea.lp^l was applied as well, but 
considering the fact that the wedge product disappears 
for k = k'. 

Eq.JSU emphasizes the wedge product kAk'. As to 
be seen below, it is closely related to the orbital angular 
momentum L, which likes for example to be dotted into 
the spin to produce spin orbit coupling L-er. It is therefore 
reasonable to divide the spinor factor into two parts, 



S — So + S'l 



(35) 



with Sq being independent of the wedge product, and with 
S'l depending only on the wedge product. Thus: 



(s 1 ,s 2 \S \s' 1 ,s' 2 ) = 



E 1 (k)+m 1 E 2 (k)+rn 2 
1ra\ 2m 2 



(36) 



J E l( k^ +r E 2[ k^ + 
k-k' 



_i k-k' | k-k' 

~T (Bi(fe)+mi)(Bi(fe')+m 1 ) + (E 2 (k)+m 2 ){E 2 (k')+m 2 ) 

, k 2 k' 2 

+ {E 1 (k)+rn 1 )(E 1 (k')+m 1 )(E 2 (k)+m 2 )(E 2 (k')+m 2 ) 

k 2 -(<TiAk)-(g 2 Ak) k' 2 -( CT iAk')-( CT2 Ak') 

+m 1 )(E 2 (k)+m 2 ) + (E 1 (k')+m 1 )(E 2 (k')+m 2 ) 

kk' + (<ri Ak)-(er 2 Ak') - 



K. —yoi/\H.)-yo 2 /\t<*) 
{E 1 (k)+m 1 ){E 2 (k)+m 2 

kk' + ((TiAk)-(cr 2 Ak') 



kk +((TiAk)-^cr 2 Ak J kk +(<ti/ 

+ (E 1 (k)+m 1 )(E 2 (k')+ m2 ) + (B 1 (fc')+mi)(B 2 (fe)+m2) 

and 



'J ri 



) S 2)j 



(si,s 2 \Si\s[,s 2 ) = 



r\_ I E x (k)+ mi E 2 (k)+m 2 



2m i 



(37) 



E 1 (k')+m 1 E 2 {k r )+7n2 
2mi 2m2 



(si,s 2 \ [ 



i(Ti -kAk ' 



i(72 -kAk 



+ 



(Ei(fc)+mi)(Ei(fe')+mi) 1 (E 2 (k)+m 2 ){E 2 (k')+m 2 ) 

-(kAk') 2 -(o-i-kAk')(<j 2 -kAk')+i(cri+ ( T2)-kAk'(k-k') 

(E 1 (k)+ mi )(E 1 (k')+m 1 )(E 2 (k)+m 2 )(E 2 (k')+m 2 ) 

i(cri+o- 2 )-(kAk') i((Ti+(T2)-(kAk') -, , 

(B 1 (fe)+mi)(B2(fc')+m 2 ) + {E t (k')+ mi ) (E 2 (k)+m 2 ) J I' 



'is S 2)- 



One should note that Si can be omitted if one is interested 
only in a spherically symmetric wave function. 



6 The Fourier approximation 

Of course, one can diagonalize the Hamiltonian directly in 
momentum space, see Ea. (|27[) . But momentum space does 
not appeal to intuition, at least not to my intuition. One 
would like to Fourier transform to configuration space, in 
order to interpret and to understand what one is doing. 
This however turns out to be impossible for the general 
case, for two reasons: (1) The Fourier transform can not 
be taken analytically, and (2) the square root in E{k) = 
\J m 2 + k 2 induces strong non-localities in the potential 
energy part of the Hamiltonian, see for example App. [Dj 
In the further analysis, I will apply therefore what I 
call the Fourier approximation 



E(k) 



(38) 



But one must be careful. The un-considered application 
of Ea. H38|) for example to the l.h.s. of Eq.J3U|) gives plain 
non-sense. 

The non-sense is avoided by substituting in Ea. H27|) 



A(k) 
B(k) 



(39) 



Similarly, one substitutes the four-momentum transfers in 
Eq.(J2J) according to 



(k- 
(k- 



k') 2 
k') 2 



{Ei(k) 
(E 2 (k) 



E^k')) 2 
E 2 (k')) 2 



(k- 
(k 



k') 2 
k') 2 



(40) 



I try to avoid the expression "non-relativistic approxima- 
tion." Relativistic or non-relativistic motion is usually as- 
sociated with the kinetic energy. But here the Fourier ap- 
proximation takes place essentially in the interaction, un- 
der the integral, and is applied, as mentioned, for the sole 
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purpose of making a Fourier transform possible. Whether 
this is a good approximation or not, whether Ea. (|39|l is 
applicable or not, can be answered only in the solution, a 
posteriori. In any case, this question must remain on the 
agenda. 

In the Fourier approximation, Eas. (|36l37[) become very 
much more simplified. Suppressing explicit reference to the 
bras and kets \s\, s 2 ) one gets simply 



S = l 



-(<riAk)-(<T 2 Ak) , k^-^iAk^-^aAk 7 ) 



4mim2 



Am 1 7712 



kk' + ^iAkj-^Ak') k . k / k . k / 

2m\m-2 4m'i 4m'i 



J_ k^ k'^ 

16 m? mi 



q tiTi-kAk 

m — — t— 1— 



io~2 -kAk 
7^— 



2(q-i+q- 2 )-kAk / 
2mi mi 



+ 



4mf 4m 
kk 7 t(<n+<r 2 )-kAk' (kAk ') 2 + (cn -kAk ') (<r 2 -kAk ') 



(41) 
(42) 



4mim2 4mim2 



(4mim2) 2 



respectively. 

With the usual hyperfine approximation, see for exam- 
ple 0, 

(o"i A b) ■ (<72 A c) =>■ § (0-10-3) (be), (43) 

with b or c being any of the vectors k, k', or k A k', these 
equations can be simplified even more. One gets 



So = l + (-o-icr 2 )£^ 

u iy * 6mim 



(k+k') 2 




hyperfine term 
kk' 



4m±m 



(44) 
(45) 



The last term must be suppressed for consistency with 
Ea. ljSgl) . Obviously, the identity 



1 



1 



1 



mi 
m 2 



m 2 
mi 



mj m 2 mim2 
has been applied here. Correspondingly one gets 

Si 



i<7\ -kAk' 


1 _ 




1 zo"2-kAk' 


1 




4m 1 


m r 


mi 


' 4m 2 




m 2 



(46) 



spin— orbit tcrn^ 



spin— orbit term 2 



-<T 1( 7 2 ( kAk Q , ^(Q-l+^j-kAk 7 kk' 

6mim 2 4m\m<2 4mim 2 4mim 2 



(4mim 2 ) 2 

L 2 — term hyperfine— L 2 — term spin-orbit Darwin term 
= Ss 0l + S' S02 + Sl + 5hfL + SsoD- (47) 

The nomenclature is more transparent in configuration 
space. The kernel of the potential energy is then 



tf ( k;kO=-^^S(k;kO 
with q = k k' and S = So + Si . 



(48) 



6.1 The hyperfine approximation 

Let me first investigate 

(0- A b) A c = (a ■ c)b - (b • c)a . (49) 

The 3-component of Eq. (|49|l is then 

[(a A b) A c] 3 = (<7i ci + cr 2 c 2 + 0303)63 
- (61 ci + 6 2 c 2 + b 3 c 3 )a 3 
= ((7i ci + a 2 c 2 )b 3 - (&1C1 + 6202)0-3 . 

If one has reasons to believe that the off-diagonal com- 
ponents of biCj cancel by symmetry considerations, i.e. 
(biCj) ~ ^Sijhc, one gets ((crAb) Ac) 3 ~ — |(73(bc). The 
other components behave correspondingly. 

Since (om A b) • (a 2 A c) = - [(01 A b) A c] ■ a 2 , Ea. ftty 
is valid approximately. One is accustomed to this substi- 
tution from the theory of hyperfine interactions JS], and 
I will refer to it as the the hyperfine approximation. As 
shown in the appendices, the hyperfine approximation is 
a rather weak assumption. Sometimes it is even exact. 



6.2 Deriving the Singlet- Triplet model 

The wedge product (01 Ak) • (a 2 Ak') is defined in Ea. (p?2*|) . 
Replacing it by 



((7! A k) • (0-2 A k') ((7K7 2 ) (kk') , 



(50) 



rather than by the hyperfine approximation in Ea. H43l) . 
one gets from Ea.l|41[l in a first step: 



So = 1 



kk' , kk' 

4m 2 ' 4m 2 

k 2 (l-g-iCT 2 )+k' 2 (l-CTia 2 )+2kk'(l+g-io-2) 
4mini2 



Adding and subtracting 2kk'(l — o\o 2 ) gives 



= l + (l-(7l(7 2 ) 



(k-k') 2 


1 kk' 


1 




4m 1 m 2 


r 4 


mf 




Q 2 4 


kk' 




2 


4mim2 


4 




mim 2 



(51) 



In previous work pQ, the matrix S was replaced by its 
diagonal elements, 

(s 1 ,s 2 \S\s' 1 ,s' 2 ) => 5 Sl;S / (5 S2;S / (si,s 2 |5|si,s 2 ) • 
Ea. (|51|) gives then, up to terms proportional to kk', 
(si,s 2 |5|si,s 2 ) 



Q 2 



Q 2 



1 — - — , for si = — So, 

4mim2 ' ' 

for Si = s 2 . 



One has thus derived the Singlet- Triplet model ^] with- 
out ad hoc procedures. The model has played an impor- 
tant role in the development of the theory, particularly its 
renormalization. 
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7 The Hamiltonian in configuration space 

The eigenvalue equation i|25|l in momentum space, 

/ ( i 3 k'( Sl , S2 | J H'(k;k')|s / 1 ,4)¥'44(k')=£;¥' sl « 2 (k), 

can be Fourier transformed to configuration space, 

Jdh\ Sl , S2 \H(r;r')\ S [,s' 2 )V s , s ^r') = EV SlS2 (r). 

The eigenvalue equation 

H& = EW, (52) 

has then the structure of the familiar Schrodinger equa- 
tion. Quite in general, working in configuration space has 
the advantage of being closer to conventional quantum 
mechanics and to phenomcnological models where our in- 
tuition comes from. 

Fourier transformations need a sign convention, 



H(r: r') = J d\ e* kr j |^ e ~* kV iJ(k; k') , 
^ SlS2 (r) = jd\ e ik V SlS2 (k). 



(53) 
(54) 



It could be a source of endless confusion that H (k; k') 
and H(r; r') are denoted by the same symbol H . But this 
should not be a problem since the meaning is usually clear 
from the context. 

Following the nomenclature of Eas. l|45l47|l the Hamil- 
tonian is decomposed into 

H = T + V + Vi , with (55) 

V = V +V M + V K + V D , (56) 
Vi = V SOl + V SQ2 + Vl + Vml + V soD . (57) 

The kinetic and the potential energy are respectively 
„ 2 



T = 



V 



2m r ' 



* „ iqr «c(g) R(q) 
q 2tt 2 g 2 



The hyperfine, kinetic, and Darwin potentials are 

6mim2 



V K = 
V B = 



V 2 
— P , 



/ V 2 
I 4roi 



m2 1 16mim2 



mi I m 2 
ni2 171 1 



(58) 
(59) 

(60) 
(61) 
(62) 



respectively. The spin-orbit, L 2 , hyperfine-L 2 , and spin- 
orbit-Darwin potentials are 



V — CTlL 



Vl = 

VhfL = 







1 d.V 


m r mi 




r dr 



1 d z V 

dr' 2 



(4mim2 ) 2 

q~l<72 

6mirri2 4mim2 



1 <fy 



t/ _ (q-i+g 2 )-L 
^ soD (4mim 2 ) a 



( 


i dy" 


V 2 - 


'l d(V 2 V)~ 


) 




4r 




r dr 





(63) 
(64) 
(65) 
), (66) 



respectively. They depend on the angular momentum op- 
erator L = r A p, see App. [UJ and vanish for spherically 
symmetric wave functions. The linear momentum opera- 
tor is denoted by p = — iV- 

Some of the above terms appear also in the analysis of 
the Dirac equation with an external vector potential, see 
for example One should emphasize that the present 
analysis is more complete since retardation effects are fully 
included. They generate the remainder of the terms in the 
above. 



8 Spherically symmetric wave functions 

It is our freedom to solve Ea. (|52(l only for wave func- 
tions with spherical symmetry, if r SlS2 (r) = ^P SlS2 (r). The 
s-states, in a way, are also the most interesting, since prac- 
tically all hadrons suitable for beams or targets have this 
symmetry. And for them one can solve the problem rigor- 
ously. 

Since there is no angular momentum, one can couple 
the quark Bj0rken-Drell spins into total spin and thus total 
angular momentum, by means of Wigner-Eckart |12j , 

\ S > S >) = E^H^Mh), (67) 



with the total spin (or total angular momentum J) is ei- 
ther S = for singlets or S = 1 for triplets. Using the 
completeness relations of the Clebsch-Gordan coefficients 

J2s,s z (l s ^ s 2 \S,S z ){S,S z \±s[± 4) = 5 S1 , S '/ S2 , S ^8) 
I can perform (for the last time) a unitary transformation 
£s,s.(k) - £ vWk)<i Sl i S2 |S,^) , (69) 
V s , Sx (r) = ^s 2 (r)(^s 2 \S,S z ), (70) 

S1.S2 

and solve the Hamiltonian equations 

H$ s ,s z {r) = E$ s ,sAr), (71) 
or H ® £ S ,S, (k) = E$s, Sz (k) , (72) 

for fixed and given S and S z rather than carrying out the 
spin summations in Ea. H52(l explicitly . In this represen- 
tation, the only spin-off-diagonal operator in Ea.l j60|l is 
diagonal, i.e. 

™-m(s + i)-3-{5 £|:j;S3ic«) 

It depends on the taste, whether one works first in con- 
figuration space and subsequently Fourier transforms to 
momentum space, or whether one solves directly Eq. l|72|l 
in momentum space. At the end one has !ps,S z O s -)i which 
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can be unitarily transformed back to the light-cone wave 
function according to 



^ A(k z (x), kj_) 



y/x{l-x) 

x </nfc 2 ]fl t (fc*(a:),kj L )|s 1 a 2 > 

si,s 2 

s,s 2 



(74) 



see also Eas. (|l(Jf) and (|21|l . Note that the light cone wave 
function is a superposition of singlet and triplet. 



9 Perspectives 

I have promised a technical paper. Sometimes technicali- 
ties are important. 

Most of the above work relates to back ground knowl- 
edge on atomic and Dirac theory and the results speak for 
themselves. 

Instead of discussing the results to some detail, I will 
highlight only some important aspects: 

• Most remarkably, the fine and hyperfine interaction 
depends on the coupling constant only through the 
potential energy V(r). 

• One faces a seemingly non relativistic potential energy 
which is directly related to the light cone wave func- 
tions with all their wonderful advantages to calculate 
cross sections, structure functions, distribution ampli- 
tudes, and the like. 

• In fact, one has the wonderful advantage of applying 
phcnomcnological approaches working with potential 
energies, and transforming their solutions to light cone 
wave functions. One can then calculate dynamic quan- 
tities like cross sections. Thus far, this was possible 
only with a lot of hand waving. 

• This opens a broad avenue of model tailoring and com- 
parison to experiment. 

• It seems one has gotten the cookie and the cake. 
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A Fourier transform of the kinetic energy 

The kinetic energy is defined in Ea. (|26() 



T(k;k') 



2m, 



-<5( 3 )(k-k') 



The Fourier transform according to Ea. 1)53(1 is 



|d¥T(r;r')!?(r') = 



(75) 



<tV lePk j dV 6 ' ^ Y T(k; kO^CrO 

k 2 ik(r-r') _Y7 2 

- , " ; •■ /A 2^W S ' (t ' ,= 2^ S ' (t, ■ 



B Fourier transform of the potential energy 

Inspection of Ea . I(45|l and Ea. ((47|l reveals the potential 
energy t/(k; k') to depend on the momenta through the 
combinations q = k — k', (kk') and (kAk'). It is therefore 
convenient to introduce sum and difference, i.e. 



k =P+f, q= k k', 
k' = P §, p=i(k + k'). 

The typical combinations are then 



(76) 



k.k' = P 2 -a_, 

k A k' = q A p , (77) 
kr - kV = p(r - r') + f(r + r') . 

The general Fourier transform Eq. (|53|l is then replaced by 
H(v;v') = J^J^e^-^e^+^H^p), (78) 
i.e. the kernel must be expressed in terms of q and p. 

B.l The central potential 

According to Ea. (|45|) and (|47fl the central potential is 

a c (q) R(q) 



U c (q;p) = U(q) 



q 2 2tt 2 
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Define the function V(r), 1. Consider the kernel of a general integral equation: 

[/ s (k;k') = /(k ) k')[/(|k-k'|) (81) 



(82) 



V(r) = d\ e l ^U{q) , (79) 

J 2. Introduce p and q according to 

and get k =p+f, q= k-k', 

, , e i P (r-r') k' = P -f, P=i(k + k')- 

V c (r;r') = Jd\ e^ r+r >£%) J^P-^- . 3 . Express U g in terms of p and q: 

= <5( 3 )(r-r')V». (80) C fl (q;p) = F(q,p) U{q) (83) 

Folding with the wave function, 4. Substitute q and p according to: 



V c ®^= Jd 3 r'V c (r;r') <P(r') = V(r)&(r) = VW , 



q => -«V y , 
p => -iVf. 

5. Get the folded Fourier transform by 



(84) 



produces the potential energy V(r) as a multiplicative fac- 
tor. V % ®W = F(-iV v ,-iVw)(V&). (85) 



B.2 General rules 

Consider first the simpler case of J7 a (q;p) = qU(q). Ac- 
cording to Ea. (|8t)|l its Fourier transform is 



This suffices to Fourier transform all functions of interest. 

C The different terms of the interaction 

According to Eas. i|45l47fl the hyperfine term is 



K(r; r') = Jd\ Jd^p ^(q f^'^y^ ■ U ^ k > = ( k " k '^ U ^ ^7 > ( 86 ) 



(2tt) 

Integrating first over p gives a Dirac delta function 



C/hf(q;p) = q 2 U(q) 



6mim2 ' 

m 2 



V a (r; r') = <S (3) (r - r') J d\ qU(q) e iqr . The kinetic term is 



M^) = £f/(9), (87) 

„2 



The q can be substituted by the partial q — iV r 

which gives K(r;r') = ^ 3 )(r - r') (-*V r V(r)). Folding fkfap) = 2> 2 TO| ,„, 

gi ves o cvan 

Ihe Darwin term is 

where Vy acts only on V(r). It is clear how to generalize TT /i,. u\ _ kjc' r7Y«1 

4.1. • 4- u-4. tt / \ I \n Tit \ lVd(.k, K J - jr^-U(q) , 

this to arbitrary powers t/ a (q; p) = (q) t/ (?)• rs 

Next, consider E/ b (q;p) = pU(q). Take the Fourier U D (q;p) = {<f - jP 2 ) , 

transform first over q, thus , r " 

J(p(r— r ) 



d3fJ P (2tt)3 7 ^ C/ ('?) eiq(r+r ' )) The typical spin orbit term is 

lW(k;k') = ^^£/( 

f/ SOl (q;p) = ^pf/(g) 



F(i(r + r')) fa> p e ' ( . P(r :;' ) . C/ SOl (k;k') = ^^^) [£ + , (89) 

J (27TJ 



m r mi 



m r mi 



The p can be substituted by the partial p =>- — iV r , 

which gives Vb(r;r') = V(r) (™iV r <5 (3) (r - r')). In the V sol ® <P = - iov J^ AV|i> (TAP) 
folding, the derivative of the delta function generates a 
derivative of the wave function <P(r), The L 2 -term is 

Vh ® IP = -iV»(VB>) , ^(k. k0 = _^^Z7(«) , (90) 

where VV acts only on !P(r). It is clear how to generalize t/L(q; p) = — ^ v ^ x2 U(q) , 

this to arbitrary powers f/&(q;p) = (p) n t/(g). mim 2 

These findings can be cast into the general rules: Vl <8> !P = ~^^7m^P" (^^) ■ 
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The hyperfine-L 2 term is 



C/hfL(k;k') 
%fL(q; p) 

VhtL ® $ 



(kAk')' 



_ (qAp) 2 ^(g) 



(91) 



It is plainly impossible to find an analytical Fourier trans- 
form of this, except when applying series expansions in 
k 2 /m 2 , i.e. 



_ (VyAVaQ 



4mi m2 

The spin-orbit Darwin term is 

r l\ _ i(o-i+o-2)-kAk' 

^oD(q;p)- i(<T1+g2) - qAp '" 2 

VsoD ® $ 



6mi7Tl2 ' 



C/ soD (k; k') = h^j-kak (k , k/) [/( j i_ 



(92) 



4mim2 W 
'/(cri+fT2)- Vv A W 

4m 17722 



finally, which completes taking the Fourier transforms. 

The potential V(r) is strictly spherically symmetric. 
With 



r ar 



(93) 



and with the usual orbital angular momentum 

L = -ir A V , (94) 
one can simplify the above equations considerably. With 

ldV- 



(Vy A V*) (F*-) = (VV A V#) = 
one can substitute in Eqs. (JH2J)— (|S3> 



(Vy A VsO 
and replace them by 



l_d 



rAVf 



m r mi 



g" i -Lip 

4mi 

i d 2 



l_d_ 

r dr 



V so D ® & 



0\0<l (Lj>) 
6mim2 4mim2 



1 eT 



(95) 
(96) 
(97) 



(4mim 2 ) 

a very suggestive form indeed 



I A. 

r dr 



(|V| - V 2 V ) V¥, (98) 



E(k) =m\ 1 + — 



k 2 

2 



Hi 



1 

2 \ m 



k 2 

2 



8 I m 2 



(100) 



but then all the beauty of the present approach gets lost. 
In the Fourier approximation, Ea. (|38|l . only the first term 
is included. 



D The non-local central potential 

One can define the kernel of a spherically symmetric po- 
tential in which the Fourier approximation has not been 
made. Restricting to the '1' in Ea. (|23|) . one gets from 
Eq.lEZl): 

m 

Zir' i LJ z 4mim2in r 
« / Mk)E 2 {k) 

(£/i(K)+mi)(£2(fc) + m2) 



E 1 {k)+E 2 {k) 



